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ABSTRACT 


This paper presents the results of efforts to implement a hybrid 
control system at the U. S. Naval Postgraduate School. The salient 
features of this effort are use of real tíme digital computor pro- 
gramming, digital determination of instantaneous state variables, and, 
using the foregoing features to demonstrate optimum control of a 
second order system. It was found that the techniques employed would 
control the second order system, and that the derivative of a variable 
could be determined with the digital equipment. 

The computer programs for accomplishing this effort are included 


in Chapter IV and the appendicíes. 
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CHAPTER 1 


INTRODUCTION 


The purpose of this work was to implement a hybrid control system 
at the U. S. Naval Postgraduate School. 

The control problem was approached assuming that a prospective 
control policy would require instantaneous knowledge of state variables 
which could not be obtained from the analog sub-system. States needed 
to implement the control policy would therefore be calculated by the 
digital computor. 

In order to achieve control and obtain the states (i.e., deriva- 
tives) it was necessary to be able to determine with congiderable 
accuracy the times required to sample and to compute. Also, appropriate 
constant time floating point arithmetic and converting procedures had 
to be developed. 

In addition, it was desired to determine what effects on state 


determination and control any noise present in the system might have. 





CHAPTER II 


EQUIPMENT 


The system under consideration is shown in Fig. 1. The digital 
computer was a Control Data Corporation 160. The memory capacity is 
4096 twelve bit binary word with a storage cycle time of 6.4 
microseconds. Operation is controlled by а program internally stored 
in sequential locations. More complete information including a 
listing of computer instructions is contained in Ref. 1, 

The plant was simulated on a PACE TR 10 Electronic Associates 
Incorporated analog computer. It is a twenty amplifier computer with 
a saturation level of + 10 volts. The power supply is regulated to 
+ 10 millivolts. Further information may be found in Ref. 2. 

The A/D and D/A converters were designed and constructed at 
the U. S. Naval Postgraduate School using Digital Equipment Corpora- 
tion modules. The principle of operation is successive approximation 
conversion. The analog capacity was O to -10 volts with a sensitivity 
of 2.4 millivolts per bít. The digital ranges 0000 to 2777 and 4000 
to 7777 correspond to analog ranges of =5 to 0 and -10 to -5 volts 


respectively. 
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FIGURE 1 





CHAPTER 11] 


DICITAL TIME SCALING 


Of the factors to be considered in the implementation of a hybrid 
control system one of the more important is the time scaling of the 
digital computer. Normally, time is not one of the criterion of digital 
computer programming, but when the computer is used as part of a control 
system, it must be time scaled to be compatable with the remainder of 
the system. Often this time scaling is referred to as computing in 
"real" time; that is, the time scale of the analog system. The follow- 
ing discussion presents some of the aspects of digital time scaling. 

All digital programs must be written so that the time required 
for the program is independent of the magnitude of any input or computed 
variables. In many cases this involves padding the shorter paths to 
equalize all path times to that of the longest path in a particular 
operation. As a note of caution, care should be exercised to avoid a 
prohibitive delay which would negate the value of the program. 

Any library routine used in a compensation program must meet this 
constant time requirement. In most cases thie was not one of the original 
requirements, meaning the routine must be rewritten. This hes been done 
for the floating point add, subtract, multiply and divide and the fixed- 
to-floating point and floating-to-fixed point conversion routines. 

These are presented in App. I as the Constant Time Floating Point 


Arithmetic Package. 





In writing constant time programs for the CDC 160 the cycle time 
was used as a basis. Each of the CDC 160 instructions requires one 
to three cycle times for execution. By writing programs so that all 
possible paths contain the same number of cycles the constant time 
requirement was satisfied. Knowing the cycle time, the program time 
was calculated. 

Although Ref. 1 states the cycle time for the CDC 160 is 6.4 
microseconds, it was found that the actual time may vary slightly 
from that value. In the case of programs which contain several 
hundred cycles the difference may become significant. Ву timing the 
below program, which contains 20 x 100 cycles, with a stopwatch the 


cycle time may be determined. 


1000 2600 icc 00 1010 2200 1de 00 
1001 2341 2341 1011 4512 4512 
1002 4213 etf 13 1012 0701 өз 01 
1003 2600 1сс 00 1013 6501 nzb Öl 
1004 7640 17640 1016 7700 hie 
1005 4211 stf 11 1015 0000 

1006 5610 aof 10 1016 0000 

1007 6501 nzb O1 


As another consequence of the constant time requirement, any 
satellite equipment used with the computer must be tested for operation 
time. To illustrate, consider the analog to digital converter, There 
is a time delay between the time an input is called amd the time the 
input is inserted into the computer. During this delay the computer 
cannot be used for other purposes. This time delay must be found so 
that it can be taken into consideration in computing total program 
time. To determine the delay, run a program of the input routine ree 


peated several times and with a stopwatch time the program. Knowing 





the total number of analog to digital conversions and the time for 
these conversions, the time per convereion can be calculated. A 
similar procedure can be used to determine the digítal to analog con- 


version time. 





CHAPTER IV 


HYBRID CONTROL POLICY 


1. Theory 


The process description is 
6 = F e(K) * D U(K) IV-1 
where F [95] and D = [2] р 
and in sampled form is 
Ө (К+ 1) = еРТө(к) + 2 8009) 10-2 
Let а = еРТ = I + FT + (FD? + 
ë = I + PT + (P) nz — 


gives @ = Lar] IV-3 


and A= [Т «Ора 
-£ BOE e = [ra [4] we 
Substituting into the process equation, we obtain 
е(К + 1) = ах (К) + U (K) IV-5 


It is now desired to get the reference position in 2 samples 


8 (1) = 48 (0) + U (0) 

Ө (2) = 490 (1) + UG) = O 

ө = & (&X (0) + U(0))+ UC) 

ә -fè r1 X (0) + [2272] 0(0) + Ба U(1) 1У-6 
Letting U, = -& Ө (а) -а) 6 (п), and IV-7 
expanding IV-7 

е (0 + 2T 8 (9 = -2 Two - Z vq) IV-6a 

9 (0) = TU(O) + т ш) (x 1⁄2) 


IV-6b 
Simultaneous solution of these equations gives 


سی وھ = ко‏ 


or generally, U(K) = ар = 2 К IV-8 
the desired control policy. This policy will force the second order 
system to its reference position in two sampling periods (2T). (The 
sampling period (T) should not be confused with that alluded to in 
the state generator discussion in App. II) 
2. Memory Required 

All lower memory, except 70-77. 

Upper memory - 1645 celle. 
3. Subroutines Required 

Floating Point Arithmetic Package. 

State Generator Program with the "Divide-by-Four” subroutine which 
it requires. 

4. Special Features 

The sampling period (T) is adjusted by manually entering the 
desired counting value in cell 2430 (ae the program is herein listed). 
1712 for T=1 second, and 730 for T=} second. 

When the sampling period T is changed, the equivalent floating 
point number must be entered into cells 50 and 51, and its square into 
cells 52 and 53. These entries give the proper gain value for the 
chosen sampling period. 


5. CDC 160 Computer Control Program 


The program to implement the control policy follows. 





Summary of Cell Allocations 


(including detail of lower memory) 


0100 - 1703 floating point arithmetic package 
1704 - 1730 Divide-by-4 

1731 - 1760 2T/3 computation 

1761 - 2303 State generator 

2304 - 2414 Control Policy including output 
Lower Memory 

01-04 floating point arguments 

06 program loop-closing jump (to 1761) 
07 subroutine return jumps 

10-1 UW & L floating point forcing fet 
12-3 TE2U € L 36/27 

14-5 ТЕ10 & L 0/T? 

16-7 21030 « 1, 3/2Т 

20 U forcing function 

21-2-3 ӨА, ӨВ, ӨС 

24-5 ӨВ-ӨА, ӨС-ӨВ, resp. 

26 (9C-0B) - (98-9А) 

30-1 (ӨВ-ОА) 0 6 1, 

32-3 OBDU & L 

34-5 (9C-9B) U € L 

36-7 @CDU & L 

40-1 (ӨС-ӨВ) = (ӨВ-ӨА) / 61. 

42-3 есори & L 

44-5 “4 TU «€ L 

46-7 A TZU € L 

50-1 TU & L 

52-3 T2U € L 

54-5 2U € L 

56-7 3U & L 

61 Divided by 4 (1704) 

62 FSB 0317 

63 FAD 0335 

64 FMU 1315 

65 FDV 1533 

66 FIX-FLO 0100 

67 FLO-FIX 0162 

Unused 


00,05,27 € 60 





E 
1732 
1759 
1734 
ШЕР 
1156 
317 
1740 


17+1 
1742 
1743 
1744 
1745 
1746 
1747 
27/50 


2751 
1752 
22 
1754 
25/55 
1756 
ПЕРО 
1760 


Continue into the main loop 


2054 
001 
2055 
1002 
2050 


4003 


2051 


Look 


0101 
0604 
4007 
7064 
2056 
4005 
2251 


hooh : 


0101 
0604 
4007 
1065 
2001 
4016 
2002 
1017 


FORM 2Т/З 


10 


of the 





2504 
2505 


* 2306 
2307 


2310 
2311 
2312 
2313 
2514 
2515 
2316 


22217) 


2320 
2521 
2322 
2323 
2324 
2325 
2326 
2327 


2330 
221 
2332 
2555 
2334 
2222 


- 2336 
Bu 


2340 
234] 


2 254 


2343 
2344 
2242 
2546 
2347 


0101 
0605 
LOOT 
2023 


7066 
2052 
4003 
2093 
4004 
0101 
обон 
4007 


7065 
2001 
LOLL 
2002 
4015 
2036 
4001 
2051 


002 
2016 
1003 
2017 
коон 
0101 
0604 
4007 


7065 
2001 
4012 
2002 
4015 
2014 
4005 
2015 


pta 


std 


144 


Jpi 
таа 


199 


зта 
pta 


adn 


std 


Jpi 
144 
std 
144 


18а 
std 
144. 


std 
144 
std 
ea! 
std 
pta 


std 


Әрі 


gd 


std 
lad 


lad 
std 
lad 


CONTROL POLICY 


ӨС to floating point 


О(2 25 


© DOT/ (21/3) 


ET 





е. 


2550 
2551. 
2352 
297) 
2554 
2555 
2556 
2557 


2560 
2361 
2362 
2565 
2364 
2565 
2566 
2567 


2570 
201 
2372 
2513 
2574 
8519 
2376 
2371 


2500 


. 2h01 


2402 
2h05 
ol ok 


2405 


2406 
2407 


2410 
2411 
2412 
2413 
2414 


Look 
0101 
0604 
4007 
7063 
2001 
4010 
2002 


4011 
0101 
обон 
4007 
7067 
1600 
ШЫ 
4020 


7500 
2401 
7304 
0021 
6103 
6002 
0020 
2200 


17242 
4211 
2200 
0117 
0701 
6501 


0701 


6507 
6202 
0000 
7006 
7717 


stà Oh 
pta 

адл ОЦ 
std 07 
¿pi 63 


194 01 


std 10 


1995025 


stû 1l 
pta 
adn Oh 
std 07 
[7 
sec 
{ЧИМ 
sta 20 


exf 00 
2401 
out Oh 
0021 
net 05 
zjf O2 
0020 


14с 00 


1112 


STILL 


laf 00 
0117 

sbn 01 
nzb 01 
laf ol 
sbn 01 


nzb 07 
pur 02 
0000 

¿pi 06 


horae 


840 to form -U 


-0 to fixed point 


complement gives 40 


output U 


jump to close main loop, (0006)=1761 
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CHAPTER V 


PROCEDURE 


The programs for generating the state variable (App. II) and 
control function (Chap. IV) were tested using the circuits shown in 
Figs. 2 and 1, respectively. 

The TR-10 biasing levels, Fig. 2, for closed loop control were 
very sensitive. A small error yielded completely erroneous results. 
This was not a problem in the state variable test. The limiter was 
provided to protect the A/D converter from inputs beyond its range 
(O to -10 volts). The linearity of the A/D and D/A was determined 
to be satisfactory before operating the system, Initial conditions 
were chosen to avoid saturating the voltage limits (e.g. U(O) = ө,/т2); 
also, saturation of the float-to-fixed poínt program will occur when- 
ever the floating point number exceeds the largest possible positive 


and negative fixed point numbers, 3777 and 4000, 
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FIGURE 2 





CHAPTER VI 


RESULTS AND DISCUSSION 


The program for generating the first derivative is shown in 
App. II. Using this program and the system of Fig. 2, the average 
peak-to-peak noise level is about 2.0 volts. Except for the noise 
effects, the derivatives obtained were reliable. Any error in 
amplitude and/or phase being within the noise envelope as shown in 
Fig. 3. 

The program for implementing the control policy is shown in 
Chap. IV. The results of testing this program with an imbedded 
step are shown in Fig. 4. If noise is relatively low at the instants 
of sampling, the control is nearly ideal, as in Fig. 4a. However, 
when noise is appreciable the result obtained is like that shown in 
Figs. 4b and 4c. Fig. 4b shows the effect of noise making the Š term 
of the control policy too large, while Fig. 4c shows the effect when 


à is made too small. 
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CHAPTER VII 


CONCLUSION AND RECOMMENDATIONS 


The control policy has been demonstrated to be effective when 
noise is not a significant factor. This is seldom the situation. The 
policy as implemented here is very sensitive to noise of the level of 
several millivolts. 

It was further demonstrated that a good first derivative could 
be digitally determined, but is very sensitive to input noise. 

It was not possible to determine the major source(s) of noise, 
but the test signal generator was highly suspect. The generator 
probably has a high frequency ripple of several millivolts in its 
output, thus giving the noisy results previously discussed. 

Improvement of the results obtained here can be obtained either 
by eliminating noise or turning to more sophisticated programing, 
such as predictor-corrector schemes, Noise elimination ís all but 
impractical except in the laboratory. More sophisticated programs 
will result in much higher "costs" than the programs presented here. 
These costs will be mainly in computor usage-time and memory space 
required. Also, while the program utilized here was almost instan- 
taneous, some of the more sophisticated ones may involve significant 


time (phase) delays for which adjustment will be necessary. 
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Further investigations along the approach used here should 
initially determine the noise characteristics o£ all equipments 
considered for use in order that excessively noisy signals are not 


introduced. 
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APPENDIX I 
CDC 160 CONSTANT TIME FLOATING POINT 


ARITHMETIC PACKAGE 


The Constant Time Floating Point Package contains routines for 
fixed point to floating point and floating point to fixed point con- 
versions and floating point add, subtract, multiply and divide. 

These subroutines have been programmed so that the same amount of time 
is required for each operation regardless of the magnitude of the in- 
put quantity. The following paragraphs contain some general informa- 
tion and notes on the package and its use. 

1, Floating Point Number Format 

The floating point number consists of eight octal digits and re- 
quires two storage cells. Ome cell (lower) contains the significant 
figures of the number. To form a positive floating point number the 
fixed point octal number is written in binary form and the binary 
point shifted to the left of the first significant (one) bit. The 
octal number of shifts required added to 2000g becomes the upper 
portion of the floating point number. The lower portion is the signifi- 
cant bits in octal form, with the binary point understood to be at the 
left. To form a negative floating point number the same procedure is 
followed except the upper and lower portions of the number are "sevens 
complemented", 

Example: 


3721.% = 011 111 010 001., 


3 


3721,5 = .111 110 100 010, x 2138 


13, 


° 76424 ж 2 


А= 1 


P мэт ip ” 





ine positive floating ncint number ig 2515 and the negative 


576% 
18 ۰ 


2. Significant Figures 


Ihis package recognizes the following configurations as zero: 


Upper 0000 0000 XXXX 
Lower 0000 XXXX 0000 


In the interest of conserving time the last significant figure has 
not been rounded off in any of the routines in this package. In the 
floating-to-fixed point comversions, if the magnitude of the number 
exceeds capacity of the computer in fixed point, the output will be 
37773 if the number is positive or 4000g if the number is negative. 
3. Cells Required 
The package requires the following cells: 

01 xX upper stowage 

02 X lower stowage 

03 Y upper stowage 

04 Y lower stowage 


07 ХХХХ return jump address 
0100 to 1703 floating point package 


4, Jump Routines 
The following routines are recommended for jumping into the 
floating point package. 


Fixed point to floating point conversion: 


0101 pta XX is the number of the cell in which the 
0605 adn 05 fixed point number is stored. 

4007 std 07 

20ХХ 144 ХХ YY is the number of the cell in which the 
JOYY jpi YY address of the fixed-to-floating point 


routine is stored, 





All other operations: 


0101 pta XX ís the number of the cell in which the 
0604 adn 04 jump address of the operation is stored, 
4007 std C7 

20ХХ ірі XX 


In each case control is returned to the main program at the address 


following the jpi instruction, 


5, Summary 
OPERATION JUMP DESCRIPTION ANSWER TIME 
(CODE) ADDRESS ІМ CYCLES 

Fixed to 0100 (A) => X 01 167 
floating point 02 
conversion, 

(FX= FL) 
Floating to 0162 Х ТА А 216 
fixed point 
conversion. 

(FL-FX) 
Floating point 7 X = Y 01 319 
subtraction. 02 

(FSB) 
Floating point 0323 X + Y 01 303 
addition. 02 

(FAD) 
Floating point 1315 X x Y 01 674 
miltiplication 02 

(FMU) 
Floatíng point 1533 x. Y 01 380 
division. 02 

(FDV) 


Ar 3 





0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 


0110 
0111 
0112 
013 
011) 
0115 
0116 
0117 


0120 
0121 
01122 
0123 
0124 
012) 
0126 
0127 


0130 
0751 
0152 
0155 
0134 
0135 
0136 
0137 


0110 
ОТЕТ: 
0142 
0143 
0144 
0145 
0146 
0147 


` 


CONSTANT TIME 


FLOATING POINT ARITHMETIC PACKAGE © 


4002 
0515 
4255 
0400 
باوج با 
4001 
2002 
4247 


6204 
2402 
цоог‏ 
6204 
0101 
ہاج ہا 
2002 
6033 


501. 


5656 


6603 
2200 
2014 
3401 
4001 


2226 
6205 
2401 
4001 
2402 
4002 
2222 
5620 


4221 
6010 
2217 
0701 
4.215 
0300 
0300 
0300 


std 
len 
stf 


stf 
514 
144 
stf 


pif 
1са 
ста 
pjf 
140 
stf 
144 
zn 


aod 
aof 
srd 
pjb 
14с 


00 


2014 


sbd 
std 


Led 
pif 
led 


ed 
std 
lar 
sbf 


stf 
Zn 
taf 
sbn 
stf 
nop 
nop 
nop 


FX-FL. Store fixed point nomber in cell 02. 
SEV COLRT 714: 


Zeroize count 2 and cell 01. 


Store fixed point number for later sim check. 


Number is negative, complement enc restore. 


Add one to count 2 to аи ос positive anc 


negative petns. 
If number is zero jump fora. 


Shift celi 02 until iirst bit is acne bit R 
For each shift add one to cell Ol and 
count 1. Cell. 02 now contains the lover 
part osi Т1Ї08Х 10 number. 

Subtract cell 21 from 2014 to form the upper 

“seating point number end 


„nd. Ly La, 
= ^ EE 
остар v QUE М. E. 


- 


restore. 


Check Sign 7۰+ +) 21 
Nuber is negative, complenent upper and 


lower parts and restore. 


Ж Sr ог delay Loops 
06:20 сє 7٤ 


FORMA COLNE. Gols 


op 


NO delay 1000500666 ٦ 
Delay loops to equalize tine. 





0150 
0151 
0152 
0153 
0154 
0155 
0156 
0157 


0160 
0161 
0162 
0165 
0164 
0165 
0166 
0167 


0170 
0171 
0172 
0175 
0174 
0175 
0176 
olr 


0200 
0201 
0202 
0203 
0208 
. 0205 
0206 
0207 


0210 


8232 
0213 
0214 
0215 
0216 
0217 


0220 
0221 
0222 


шаг» 


0224 
0225 
0226 
0227 


6506 
7007 
0415 
0300 
0300 
6615 
0000 
0000 


0000 
0000 


2001: 


4260 
6305 
ook 
0701 
6501 


6005 
2401 
4001 
2402 
4002 
2002 
6051 
2001 


5600 
2000 
6046 
6346 
4240 
0714 
1237 
6202 


6302 
614] 
ОҢОО 
4255 
2650 
4227 
4402 
4627 


5626 
5623 
6504 
2222 
6006 
0404 
0701 
6501 


nzb 
Jpi 
lán 
пор 
пор 


pjb 


rda 


sbn 
nzb 


06 
07 
15 


15 


Return jump to main program. 


Form count for fixed point number 


egual to zero. 


Sign check stowage. 
Count 1 stowage. Count 1 


Count 2 stowage. 
Count stowage. 
Load upper part of f 


PL-FX ° 


= 13 - No. of shifts. 
J А 


loating point 


number and store as sign check. 


Number is positive. Time delay to equalize 


positive and negative paths. 


Number is negative. Complement upper and 


lower parts and restore. 


! 


If number is zero jump to 0247. 


Subtract 2000to find exponent: If positive 
store on shift. | 
. Exponent js Zero. Jumps to 0250" 


Exponent is negative. Jump,to: 0251. 


Subtract 14 from shift end stow result 
in count. Count = No. delay loops. 
If count is positive ,number will exceed 


capacity. Jump to 0252. 


Zeróize mask andá form complement of shift. 


Shift cell 02 shift number of times. At 


same time form mask. 


Delay loops +o equalize time. 


А-5 


T 





0230 
0231 
0232 
0255 
0234 
0255 
0256 
БЕРІ 


0240 
0241 
0252 
0243 
0244 
0245 
0216 
0247 


0250 
0251 
0252 
0253 
0254 
0255 
0256 
0257 


0260 
0261 
0262 
0263 
0264 
. 0265 
0266 
0267 


0210 
ger 
0272 
0273 
027% 
0275 
0276 
0277 


0500 
0501 
0502 
0505 
0504 
0305 
0306 
0307 


212215 
6504 
2002 
1213 
1002 


2206 


6505 
2002 


7007 
2402 
1007 
0000 
0000 
0000 
0000 
6004 


6012 
6342 
6132 
2200 
0157 
0701 
6501. 
0500 


0h00 
1007 
2200 
0155 
0701 
6501 
0500 
2324 


6303 
0۸01 
1007 
0501 
1007 
2200 
0134 
0701 


6501 
0300 
0100 
7007 
2200 
0126 
0701 
6501 


nof 15 
nzb Oh 
lad 02 
TOR 15 


147 06 


nor 95 
lad 02 


ope Oy 
led 02 
jpi 07 


zjf Oh 


2J l2 
njí 24 
па? 52 
14с 00 
0137 

sbn Ol 
nzb 01 
nop 


ldn OO 
jpi 07 
lac 00 
0135 

son 01 
nzo Ol 
пор 

ldb 24 


80۲ 
ldn Ol 
PLO 
ien Ол 


jpi O7. 


ide 3007 
0134 
son 01 


nzb Ol 
nop 

ldn OO 
Jpi 07 
ide 00 
0126 

sbn Ol 
nzo OL 


Mask cell 02 and restore., Cell 02 contains 
the fixed point number.’ | 


Check siem of original number. If negative ' 
complement and load in‘A register. If ' 
positive Load in A register. Re- 


turn jump to main program. 


Check stowage. 
Shift stowage. 
Count stowage. 
Mask stowage. 
JUmp to0253. 


Jump to 0262. 

Jump tp 0275. 

Jump to 0304. 

Floating point number is zero. Load zero 
in A register, delay and return jump | 
to main program. 


Floating point number is larger than 0.5 and 
less than 1. Load 0001 in A register, 
delay, and return jump to main program. 


Floating point number is greater than zero 
but less than 0.5. Load zero in A re- 
gister, delay and return jump to main , 


program. 


Floating point number exceeds capacity of 
computer. If number is positive load 
3777 in A register, if negative load 
4000, delay, and return jump to main 


A-6 





0310 
071 
0512 
0515 
0518 
0515 
0316 
Sor 


0320 
0321 
0322 
0525 
0324 
0325 
0326 
0527 


0550 
0551 
0552 
0555 
0334 
0975 
0336 
0551 


0540 
0511 
0342 
0343 
0344 
- 0345 
05146 
0347 


0350 
0221 
0352 
0225 
0354 
0222 
0356 
0521 


0360 
0361 
0362 
0363 
0364 
0365 
0366 
0367 


2345 
6304 


145 45 
njf o4 


2200 0 


ST 
7007 
2702 


7007 


2005 


6104 
4005 
0500 
6005 
2403 
4005 
2004 
6104 


LOOL 
0300 
6003 
arok 
YOO! 
2001 
6017 
6555 


2002 
6016 
2003 
6042 
6305 
2004 
6041 
7101 


1077 
2004 
6035 
7101 
0653 
2001 
2001 
2003 


6105 
oho2 
OTOL 
6501 
6010 
4001 
2004 
6006 


O 
7 
leb 02 
JPO 
1да 03 


nzf oh 
std O3 
nop 

95 و 
lcd 03‏ 
Std 05‏ 


144 ol, 
nzf Ooh . 


std. oh 


nop 

2105 
1са он 
sta ol 
144. 01 
20114 
Jl» 


ldd 02 
21210 
144. 05 
zjf 42 
30000 02 
1da Ok 
211 41 
(21 


1077 
laa oh 
22855 
3281 
0653 

144. 01 
144 01 
аа 5 


Диаш o 


500 01 
11203 07 
21121: 
sta OF 
jaa oh 
zjf O6 


FSB. 


FAD. 


program. 


j! á: 


Zero check Y upper and lower. If 


either are zero jump to.FAD. If both 
are non-zero, complement, restore ала 
jump to FAD. 


Sign and zero check X. 


X is positive and non- zero. Sign and 


bo ٦ 


Time delay. 


Zero check Y. 


X and Y are positive and non-zero. Jump 


X is positive end Y is negative. Both 
are non-zero. Jump to 0653. 


NOTE: X is zero 70770355 
0 0 7 


Zero check Y upper. 


Time delay. 


Y upper non-zero. Store 'in cell 01. 


Zero check Y lower. 


Dt 





0510 
0511 
ШЕТ 2 
9212 
0374 
0575 
0370 
0577 


0400 
0401 
0402 
0405 
ОЦОЦ 
0405 
0405 
0507 


0h30 
ET 
0h12 
0413 
0414 
0415 
0:15 
T 


"00 
0h21 
0422 
0h23 
Bo 
0425 
0426 
0627 


0430 
0431 
0432 
0433 
0434 
0435 
0436 
0437 


OLLO 
он 
01442 
0443 
оны; 
0145 
0446 
ош 


4002 
2200 
0215 
0300 
6106 
ОВОО 
1001 
4002 


2200 
0212 
0701 
6501 
ТООТ 


20053 


2003 
2200 


0216 
6507 
2002 
64-34 
2003 
6h10 
6205 
200h 


6411 
7101 
1105 
200} 
6415 
2401 
3405 
4275 


6077 
6277 
0614 
6056 
6356 
1600 
TTT 
4266 


ОҢОО 
1,265 
hho2 
1,663 
5662 
5660 
6504 
2402 


std 02 
14с 00 
0215 
сор 
nzf 06 


icai 00 
чёт ЖОГ 


sti 02 


14с 00 
0212 

son Ol 
nzo OL 
¿pi 07 
13d 093 


+01387 


14с 00 
0210 

EZS OF 
lad: 02 
zjo 34 


194 05 


21010 
pjf 05 
iaa oh 


28511 
ol 
LIOS 

läd Ok 
220 15 
Шат 
504 03 
stf 19 


211411 
DO 
adn 14 
zjf 56 
njf 56 
scc OO 
шин 

stf 66 


140 00 
stf 65 
sra 02 
srí 63 
aof 62 
aof 60 
nzb 04 
124 02 


Zero check X lower. NOTE: X is negative 


| | 


Y lower is non-zero. Store in cell 02. E 
Answer is now in cells Ol and (02. Delay a 
end return jump to main progrem. 


Answer 18 zero. Load zero in 'gells Ol and | 79 
02, delay and return jugo main o Ж 
, program. eT 


! 
"" "0 
5 » Ч 


t poset ' 
ти 


Y is zero. Delay zo return jump to main : ۱ T 
, program. E 


for 0h12 to Oh2h. 
Zero end sign check Y upper. 


Y is negative. Zero check Y lower. 


X and Y are negative. Junp to 1105. 
Y is positive. Zero check Y lower. 
Subtract upper parts and store result 


in count 1. NOTE: 0425, to 0652 is 
X negative and Y positive' routine. 


Exponents equal. Jump to 0527. | | 
X exponent larger. Jump to 0530. i 
Add 14 to count 1 to form shift. 

Y is much larger than X. Jump to 0511. Е: 


Complement and restore shift. 


Zeroize mask. 


Shift X lower to equalize exponents. 
Form mask. No.: of shifts = shift. 


Mask and restore X lower. 


A-8 





0450 
0451 
0452 
0453 
emn 
0h55 
0456 
0457 


0460 
0161 
0462 
0465 
01464 
0465 
0466 
0467 


0470 
0471 
0472 
0473 
ہا ہان 
0475 
0476 
0477 


0500 
0501 
0502 
0503 
0504 
0505 
0506 
0507 


0510 
0217 
0512 
0515 
0514 
0515 
0516 
0517 


0520 
0521 
0522 
0523 
0524 
0525 
0526 
0527 


1256 
4002 
0513 
5255 
2004 
3402 
1002 
6206 


2003 
+001 
0101 
0701 
6010 
5641 
4402 


6602 


2636 
2115 
2002 
4001 
2250 
6006 
oko} 
0701 


6501 
5623 
6504. 
2223 
6004 
4220 


5620 


6502 


7007 
0300 
2003 
5001 
2004 
4002 
2200 
0204 


0300 
0701 
6501 
7007 
0000 
0000 
0000 
6055 


1р? 56 
std 02 
Ten 12 
stf 53 
144 ہن‎ 
504 02 


pif 06 


пасва 
ста 01. 
län 01. 


sbn. Ol 
21210 
aof 41 
514 02 
25.02 


lef 36 
sbn 13 


ада 05: 


std 01 
Лат 50 
zjf 06 
län OL 
sbn 01 


120201 
aof 23 
nzb Ol 
laf 23 
21122015 
stf 20 
aof 20 
nzb 02 


JEN 
nop 


1290805... 


std O1 
laa Oh 
std 02 
lác 00 
0204 


nop 

sbn 01 
02 ٦ 
ро 


21129) 


' 
1 


Load -13 in count 2. 


Subtract X and Y lower and store result 
on cell 02. 


First bit ís a one bit. No shift is 
necessary. Store exponent in cell 
Ol. Time delay. 


Shift cell 02 until first bit is a one 
bit. Modify count 2 for each shift. 


Form exponent and store in cell Ol. 


zero check count 1. 


Count 1 delay loops. 


Zero check count 2. 


Count 2 dealy loops. 


Return jump to main program. 

Y is much greater than X. Load Y in 
cells Ol and O2, dealy and re- 
turn jump to main program. 


Count 1 stowage. 

Shift stowage. 

Mask and count 2 stowage. 
Exponents equal. Jump to 0604. 


А-9 


” 


ОТТ 
4304 
6244 
0400 
4306 
4404 
710 
و‎ 


9713 
6504 
2004 
1515 
4004 
0515 
4320 
2402 


3404 


1002 
6205 
0402 
0701 
6501 
6007 
Dio 


huge 
6602 
2334 
0613 
5001 
2402 
4002 
2743 


0601 
4345 
2200 
6531 
+7 
60 
2200 
0207 


0300 
0701 
6501 
ТООТ 
0515 
Laks 
2502 
3404 


sbn 14 


stb ol 
pjf 44 


stb 06 


sra O4 
srb 10 
aob 11 


aob 15 
nzb Oh 
144 ОЦ 
1ро 15 


СА Э 


stb 20. 
: .. Subtract X and Y lower and store re- 


led. 


sbd oh 


pjf O2 
làn 
sbn Ol 
20201 
21501 
aob 31 


srá 
pjb 
185 34 
adn 15 
гай 01 
led 02 
sta 02 
leb 43 


ада Ol 
stb 45 
14с 00 
6551 


stb 47. 


nzb 50 
lde 00 


2221 


nop . 
sbn Ol 
nzb Ol 
JD SO 
۔‎ 75 
stf 45 
160. 
sba Ol 


X exponent larger. Form shift. 


X is much larger than Y. Jump to 0576. 
Zeroize mask. 


| 
Shift Y lower to equalize نہ‎ S 
Form mask. No. of shifts = shift. 


Mask and restore Y lower. 


Load -13 in count 2. 


sult in cell 02. 


Time delay. 


Snift cell 02 atil TIrTSL DIL 1S 8 
one bit. Modify count 2 for each 


shift. 
Form exponent and store in се11 01. 


Answer is negative. Complement and 


restore cell 02. 
Modify count 1 to equalize delay. 


Modify program for proper jump. 
Jump back to counts 1 and 2 dealy loops. 


X much larger than Y delay. 


X and Y exponents equal. Form count 3. 


Subtract X and Y lower and store їп 
cell 02. 


' A-10 





0610 
0611 
0612 
0615 
0614 
0615 
0616 
0617 


0620 
0621 
0622 
0623 
062% 
0625 
0626 
0627 


0630 
0631 
0632 
0633 
0634 
0635 
. 0636 
0637 


0640 
0643 
обща 
0643 
общи 
0645 
0646 
0647 


0650 
0651 
0652 
0653 
0654 
0655 
0656 
0657 


0660 
0661 
0662 
0663 
0665 
0665 
0666 
0667 


цоог 
6105 
4001 
2200 
0201 
6134 
6314 
5655 


huge 
6602 
2402 
4002 
2226 
0615 
5001. 
0500 


OOO 
6013 
2402 
5002 
5616 
1502 
6602 
2615 


0715 
3003 
001 
0500 
5606 
4.205 
6502 
2200 


ОТТ Т 
6551 
0000 
2405 
3401 
14267 
6071 
6271 


0614 
6054 
6354 


1600 


0 
4260 


0h00 
1257 


std 02 
nzf 05 
std Ol 
ldc OO 
0201 

nzf 3h 
1 
aof 33 


sra 02 


pjb 
led 
std 
10% 


02 
02 
26 
15 
гад 01 


пор 


lán 
کڑے‎ 
Lead 
std 
aof 
srd 
Роб 
Lef 


00 
15 
02 
02 
16 
02 


13 
ILS 
05 
01 


06 
05 


00 


sbn 
ада 
std 
nop 
ао? 
stf 
nzb 
14с 


OTI 
NZD SL 


led 05 
съд 01 
SUE TON 
Zt 
BIE IL 


adn 1k 
zjf 5h 
njf 54 
scc OO 
ee: 

stf 60 
ldn OO 
stf 5T 


A Gnda Y are equal. Answer is с: Junp 
to delay routine. 


| ; 
X is greater than Y. Left shift cell 02 ' 


until first bit is a one bit. Modify 
count 2 for each shift. 


Form exponent and store in cell Ol. 


Time delay. 


X is less than Y. Left shift cell 02 
unail first pit is a 6 ٠٦ 
Modify count 3 for each shift. 


Form exponent and store in cell Ol. 


Count 3 delay loops. 
Time delay. 


Count 4 stowage. 

Subtract upper parts and store result in 
count I. NOPE Х 16 positive and Y is 
for 0653 to 1076. 

Exponents equal. Jump to oThT. - 

Y larger than X. Jump to OT5O. 


Ada 14 to count 1 to form shift. 


9 


X is much larger than Y. Jump'to 0735. 


Complement and restore shift. 


Zeroize mask. 


A-11 





0670 
0671 
0672 
0675 
0674 
0675 
0676 
0677 


0700 
0701 
0702 
0703 
0704 
0705 
0706 
0707 


0710 
DT 
OT12 
0713 
0714 
0715 
0716 
ОИТ 


0720 
0721 
0722 
0723 
0724 
0725 
0726 
0727 


0730 
0751 
0732 
0133 
0734 
2122 
0736 
ТІРІ 


0150 
0741 
o742 


0743 


0744 
0745 
0746 
0747 


41404 
4655 
5654 
5652 
6504 
2404 
1250 
Look 


0512 
4245 
2002 
3404 
4002 
6205 
0h02 
0701 


6501 
6007 
5634 
loo 
6602 
2631 
0713 
5001 


222) 
6006 
ohoh 
0701 
6501 
5617 
6504 
2217 


6004 
4214 
5614 
6502 
7007 
0300 
2200 
0207 


0300 
O701 
6501 
7007 
0000 
0000 
0000 
6062 


srd O 
srf 55 
aof 54 
301 92 
nzb Ol 
led O4 
lpf 50 
sta 04 


37 
stí کہا‎ 
14а 02 


съд .. 


std 02 
pjf 05 
ldn 02 
Sbn Ol 


nzb Ol 
2111 
aof 34 
srd 02 
pjb 02 
E 707 
500 13 
rad ol 


laf 2h 
zdf 06 
làn Oh 
sbn Ol 
nzb Ol 
aof 17 
nzb Oh 


tarir 


zjf OL 


COMIS 
aof 1h 


nzb 02 
ӛрі 07 
пор 
lde 00 
0207 


пор, 

sbn 01 
nzb 01 
ops Ol 


zjf 62 


Shift Y lower to equalize exponents. 
Form mask. No. shifts = Shift. 


Mask and restore Y lower. 


Load) =15 in count 2. 


Subtract lower parts and store in 
cell 02. 


First bit is a one bit: No Shift is 
needed. Time delay. ° 


Left shift cell 02 until first bit is 
a one bit. Modify count 2 for each 
Shift: | 

Form exponent and store in cell Ol. 


Zero check count 1. 


Count 1 delay loops. 


Zero check count 2. 
Count 2.delay loops. 


Return jump to main progran. 
X is much larger than Y. Delay and 
return jump to main program. 


) 


Count 1 storage. 

Shift Storage. 

Mask and count 2 storage. 
"xponents equal. Jump to 1017. 


A-12 | 





0750 
0751 
0752 
9158 
0754 
(DS 
0756 
2121 


0160 
0761 
0762 
0763 


102% 


1025 
1026 
1027 


0714 
4304 
6245 
0h00 
4306 
4402 
4710 
аг 


ВІ 
6504 
2002 
1315 
4002 
0513 
4320 
2401 


3402 
1,002 
6205 
2003 
4001 


6110 
21 


L402 
6602 
2334 
0613 
5005 
001 
2402 
1002 


2111 
0601 
L346 
0300 
6002 
6575 
6476 
2003 


4001 
2004 
4002 
2200 
0201 
0300 
0701 
6501 


sbn : 


stb 
Dis 


an 


stb 
srd 
srb 
ао 


aob 
nzb 
144 
1pb 
std 
1сп 
stb 
Ted 


ора 
std 
par 
144 
std 
144 
1121 
aob 


srd 
pàb 
140 
adn 
add 
std 
led 
std 


180 
adn 
stb 
nop 
211 
nzb 
200 
144 


514. 
144 
std 
LAC 


0201 


nop 
sbn 
nzb 


02 
02 


05 


02 
01 
05 
10 


nw 


02 
02 
34 
15 


05 


01 
02 
оа 


ul 
01 
56 
02 


13 
76 


02 


01 
oh 
02 
00 


01 
01 


y much larger than X. Jump to 1017. 
Zeroize mask. 


Shift X lower to equalize exponents. 
Form mask. No. shifts = Shift. 


Mask and restore X lover. 


Loa s 1 ٦ 

Subtract.lower varts and store result 
in cell U2. 

КУТИ + 9ٰ  ) SEIS 
necessary. Store exponent in cell 


О 
Time delay. 


Left shift cell 02 until first bit is 


a one bit. Modify count 2 for each 
shift. 
Form exponent end store in cell 01. 


Complement and restore cell 02. 


Modify count 1 to equalize time delay. 


i 


Jump to count 1 e. 2 delay loops. 


Y is much larger than X. Store Y in 


cells Ol and 02, delay and return 
jump to main program. 





1030 
1031 
1032 
1009 
1034 
1035 
1036 
1037 


1040 
1041 
1042 
1043 
1044 
1045 
1046 
1047 


1050 
TOSI 
1922 
1022 


1054 


1095 
1056 
1057 


1060 
1061 
1062 
1063 
1064 
1065 
1066 
1067 


1070 
آ0 
1072 
1975 
1074 
1075 
1076 
1077 


1100 
1101 
1102 
1103 
1104 
1105 
1106 
1107 


7007 
0515 
yall 
2hoh 
5402 
4002 
6105 
4001 


2200 
0201 
6514 
6213 
2402 
4002 
5630 
4402 


6602 
2625 
OTIS 
5001 
0300 
6212 
5620 
4402 


6602 
2402 
4002 
2213 
0615 
3003 
001 
2207 


0601 
4205 
6503 
2200 
0116 
6547 


0000 


0411 


0701 
6501 
0400 
275 
6013 
2401 
1001 
2402 


Jpi 
len 
stf 
led 
504 
std 
nzi 
stà 


lde 


07 
15 
n 
oh 
02 
02 
02 
01 


00 


0201 


nzb 
Pj 
led 
std 
ао? 
514 


pjb 
Jer 
sbn 
rad 


. nop 


3 
aof 
SIC. 


pjo 
J ed 
Sid 
161 


ада 


std. 


141 


stf 
nzb 
lac 


14 
15 
02 
02 
30 


Ver 


02 
22 


01 
ше 
02 
02 


02 
12 


05 
01 
07 


01 
05 
02 
00 


0116 


nzb 
làn 


sbn 
nzb 


stf 
211 
dca 
514 
led 


47 
цай 


Load -15 in count 3. € 

Subtract lower parts and store result 
in cell Oe. 

X and Y are equal. Answer is zero. Jump 
to delay routine. 

Х 15 larger then Y. Shift complement 


ОР пъ 1Г50 51151598 one. bit. 
Modify count 3 for each shift. 


Form exponent end store in cell 01. 


Y is Jerser then X. Shift 2211502 
Опра 5750 boit 8 8707 


Modify count 3 for each shift. 
Complement and restore cell 02. 


Form exponent and store in cell 01. 


4 


Count 5 delay loops. 


Time delay. 


ғ 


Count 3 stowage. 
Х and Y are positive. Delay and store 


sign check. 


X and Y are negative. Complement and 
restore. Store sign check. 
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1178 
ae 
1112 
1215 
uU 
RE 
1116 
E13 


1120 


1122 
1125 
1124 
1125 
1126 
1127 


1130 
11:51 
1192 
1155 
1134 
1155 
1136 
1121 


1110 
By 
nuo 
1143 
1144 
1145 
1116 
1147 


1150 
1151 
[152 
1153 
1154 
Ш 55 
1156 
Jr 


1160 
1161 
1162 
1165 
116h 
1165 
1166 
1151 


4002 
2405 
4003 
2405 
Look 
0401 
4.262 
2003 


2401 
4260 
6064 
6264 
0614 
6057 
6557 
1600 


T 
4251 
ОҢОО 
250 
Lok 
1646 
5645 
5643 


6504 
2004 
1241 
Look 
2002 
1200 
5777 
3004 


6310 
1600 


40000 


1002 
0403 
0701 
6501 
6011 


0 
0110 
0105 
ши 
2200 
2000 
1002 
5401 


std 02 
120.03 


stä 03 ' 


Led OL 
std OL 
Jan OI 
stf' 62 
1da 05 


504 01 
stf 60 
zjf 6h 
pjf 6h 
adn 14 
212221 
пат 
see 00 


TTT 

stf 51 
lan OO 
stf 50 
5га Oh 
srf 46 
aof 45 
aof 43 


nzb OL 


lád oh 


lpf 41 
stû Ok 
144 02 
16C UG 
2721 

add oh 


Nye 10 


scc OO 
4000 
sta 02 
140 05 
sbn 01 
nzb Ol 
zin 


ls 6 
TS 5 
ls 2 
1ро 15 
айс ОО 
2000 
std 02 
aod Ol 


Subtract exponents and store result in 
count. 

Exponents are equal. Jump to 1206. 

Y is latger than X. imp to 1207. 

Add 14 to count to form shift. 

X is much larger than Y. Jump to 120%. 

Complement and restore shift. 


Zeroize mask. 


Shift Y lower to equalize exponents. 
Form mask. No. of shifts = Shift. 


Mask Y lower end restore. 


Mask first bit of X lower and add Y 


lower | 
No overflow. Replace first one bit. 


Time delay. 


There is overflow. Right shift one 
place end replace first one bit. 
Store in cell 02. 


Modify exponent. 
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0070 
TJ 
1172 
Dus 
1174 
1285 
1116 
ПТ 


1200 
1201 
1202 
2205 
1204 
1205 
1206 
1207 


1210 


1 12 
1215 
1214 
215 
1216 
ШӘЛІ 


1220 
pec 
шега 
poc» 
1224 
mee) 
1226 
tee] 


1230 
1231 
1232 
1255 
1234 
1255 
1236 
07 


1240 
1241 
1242 
1243 
1244 
1245 
1246 
1247 


SOM 
6006 
040% 
0701 
6501 


6504 
бош. 


0000 
0000 
0000 
0000 
0501 
6555 
6066 
0714 


14306 
6254 
OLOO 
4310 
402 
4712 
p 
5115 


6504 
2002 
9 
1002 
2005 
5001 
2525 
0701. 


6103. 


0500 
6005 
1600 
0۷ 
4334 
0300 
200% 


1241 
3002 
6572 
0425 
0500 
0701 
6501 
2547 


Ta Il 
zjf 06 
ldn O4 
sbn 01 


nzo OL 


aof Oh 
nzb OL 
zjf 4b 


TER Ol 
nji 29 
zjf 66 


sbn 14 


stb 06 
pjf 5h 
ldn OO 
Stb 10 
srd 02 
апр ла 
aob 15 
aob 15 


nzo Оһ 
ша O2 
o 
std 02 
по 05 
std Ol 
lab 25 
sbn Ol 


120202 
nop 
ZI O, 


sce OO 


TTT 
stb 34 
nop 
laa Ooh 


lpf hi 
ada 02 
120212 
їп 25 
пор 

sbn Ol 
20:01 
lav +47 


Zero check count. 


Count delay loops. " 


Jump to delay routine. 
Check stowage. 
Count stowage. 
Shift stowage. 


Mask stowage. 


X is much larger than Y. Jump to 1262. 


Exponents equal. Jump to 1274. 
Y is larger than X. Form shift. 


- 


Y is much larser then X. Jump to 1007 
Zeroize mask. ща 


Shift K lower to equalize exponents. 
rora mask. NOs о йг ٣ 
Mask enc restore X lower. 


Store exponent r cell Ol. 
Modify count to equalize delay. 
Zero check count. 

Time delay. 


Complement and restore count. 


Mask first bit of Y lower and add 
X lower. 


Jump back to over flow check routine. 


Time delay. 


Check Sign of answer. 





1250 
ив 
12,2 
1255 
1254 
1229 
1256 
1227 


1260 
1261 
1262 
1263 
1264 
1265 
1266 
1267 


1270 
1271 
оү? 
0023 
1274 
25 
1276 
teri 


1300 
1301 
1302 
1303 
1304 
1305 
1306 
1307 


й 10 
۔‎ 1 
1512 
12212 
1314 
12112 
1516 
7 


1320 
кле 
1522 
1525 
132% 
1325 
1326 
1327 


6105 
OLOL 
0701 
6501 
6005 
2401 
4001 
2402 


4002 
7007 
2200 
0164 
6517 
2003 
4001 
2001 


4002 
2200 
Оо 
6521 
2002 
0111 
0110 
0105 


1200 
TIT 
4002 
2004 
OF 
0110 
0103 
1306 


5002 
5401 
2200 
0152 
6550 
ОҢОО 
4261 
4261 


1.261 
4261 
2-01 
6307 
4001 
2402 
4002 
0401 


nze 05 
lan ہاو‎ 
sbn Ol 


“zb 01 


237 05 
led 01 
std Ol 
led 02 


614 02 
Jpi OT 
14с 00 
0164 

nzb I/ 
1dd 03 
sta 01 
144. 01 


514 02 
14с 00 
0151 
пасе 
144 02 
16 6 
ШЕ 72 
15 2 


172.09 
aT 
sta 02 
таа 04 
15 6 
15 5 
15 2 
1ро 06 


rad 02 
aod Ol 
làc OO 
0152 

nzb 50 
ldn OO 
81101 
stf 61 


stf 61 
St O1 
еби 


137 07 


- 81 
led 02 
std 02 
140: 01 


Answer is negative. Complement end 


Answer is positive. Time delay and ян 
return jump to main program. 


restore cells O2 and Ol. Return 
jump to main program. 


X.is much larger than Y. Jump to d 
delay routine. 


Y. is much larger than X. Load Y in 
celis Ol and 02 and jump to 
delay routine. 


Éxponents are equal. Right shift 
X and Y lower one place and 
add. Store result in cell 02. 


Modify exponent. 
Junp to delay routine. 


zeroize temp. storage BEGIN FMU цөн 


, 


load first argument (upper) 
test for neg or positive 
store positive arg 

compl and store lower 


remember first number neg 
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1330 
1551 
1332 
1882 
1334 
1555 
1336 
155[ 


1340 
1341 
1342 
1343 
134) 
1345 
1346 
1347 


42550 
1551. 
1552 
35 
1354 
1222 
1356 
1521 


1360 
1561 
1562 
1563 
1364 
1365 
1366 
1367 


1370 
Ша Та 
M 
1515 
1374 
ШЕ» 
1576 
oir 


1400 
1401 
1402 
1403 
110% 
1405 
1406 
1407 


1251 
6206 
2200 
0005 
0701 
6501 
0300 
2002 


6106 * 


2200 
QD 
0701 
6501 
7007 
2405 
6507 


4003 
ہان باج 
lool‏ 
0501 
)222 
6206 
2200 
0005 


0701 
6501 
0300 
200} 
6111 
001 
4002 
2200 


0462 
0701 
6501 
0300 
ТООТ 
0518 
6505 
0000 


0000 
0000 
0000 
501 
4705 
2505 
6207 
1200 


stf 
pjf 
ide 


0005 


sbn 
nzb 
nop 
1да 


nzf 
ldc 


0477 


sbn 
nzb 
Jpi 
led 
ТІНЕ 


std 
164 


len 
raf 
рјі 
1ас 


0005 


sbn 


nzo 


nop 
194 
128 
sta 
std 
1de 


sbn 
Iz 


пор. 


ӛрі 
len 
Dos 


Stb 
srb 
1ар 
pjÍ 
їрс 


21 
06 
00 


Ol 
01 


02 


06 
00 


01 
01 
07 
07 


05 


ou 


01 
22 
06 
00 


01 


-0162 


01 
01 


07 
05 


orig number pos - delay 


zero test first arg 


repeat tests on second arg | 


jump over temp storage 


store bit count (-14) | 
shift lower ans left one 

load upper ans 

check need for carry to upper 
remove left bit from lower 
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na 


1410 
1431 
1412 
1413 
1414 
1415 
1416 
397 


1420 
1h21 
1422 
1125 
1421 
1425 
1426 
1427 


1430 
1431 
1432 
1433 
1434 
1435 
1436 
1437 


1440 
144) 
1440 
1145 
nnn 
1445 
1446 
1447 


1450 
1451 
1452 
1453 
1154 


1355 


1456 
1457 


1460 
1461 
112 
1463 
1464 
1465 
1466 
ЕНОТ 


ЭГ 
4311 
21085 
0401 
6205 
2200 
0004 
0701 


6501 
721 
2004 
6220 
2002 
5و 
1402 
0201 


6006 
2 
0501 
9232 
Оон: 
6212 
2200 
0003 


0701 
6501 
6205 
2200 
0010 
0701 
6501 
rol 


5746 
6745 
آ2 
4203 
6203 
6302 
0000 
2003 


3600 
2000 
5001 
6346 
2365 
6107 
0501 
5001 


үү 
stb 11 


aob 15 : 


làn 01 
pif O2 
1ас 00 
000% 

son 01 


nzo ol 
Sro 21 
18а он 
pjf 20 
laa 02 
таш ер 
15а 02 
lpn Ol 


231 06 
aob 32 
Шеше од 
rab 33 
ldn O1 
pp UE 
Ilac ОО 
0005 


son 01 
nzo Ol 
pjf 05 
JC O0 
0010 

sbn Ol 
nzo 0l 
sra Oh 


aob 56 
пір 45 
1ab 51 
SLT U3 
Die 75 
002 


144 03 


50с 00 
2000 . 
red Ol 
ajf YO 
lap 65 
121501 
Ten Ul 


Тай 01 


place into lower Р "E 
ехесисе Caray Dp | Ё 
set up and jump beck ' A 2 
delay 


begin next bit multiplication 
if next bit zero, no multiplication ٣ 


add to lower ans E 


check for occurrance of end around carry ú 


Jump if no corr Tor ЕАС 
corr for carry 


delay 


delay 


| 
| 
1 


return for next bit mult; if no jump, coef mult von 
restore sign record for later ref 


add exponents 


adjusts: for zero coef 


reduce exp by one 





1470 
1471 
1472 
1473 
1474 
1475 
1476 
1477 


1500 
1501 
1502 
1503 
1504 
1505 
1506 
1507 


1510 
0 
1512 
515 
1514 
2519 
1516 
1517 


1520 
522, 
Ш2 22 
1523 
1524 
1925 
1526 
1527 


1550 
1551 
1552 
1553 
1534 
209 
1556 
od 


1540 
1541 
1542 
1543 
1544 
1545 
1516 
1947 


6341 
0500 
4002 
6307 
4002 
2200 
0002 
0701 


6501 
0300 
2002 
6305 
Luge 
0501 
5001 
6205 


2200 
0003 
0701 
6501 
2336 
6006 
2h01 
4001 


2492 
4002 
6206 
2200 
0004 
0701 
6501 
0300 


7007 
0501 
ТӨЛІ, 
OLOO 
4261 
4261 
2403 
6307 


4003 
ہا ہاج 
hook.‏ 
0101 
4252 
6206 
2200 
0005 


njf 41 


dam OO ` 


sta O2 
mjf 07 
std 02 
80 
0002 

500 01 


nzb 01 
nop 

laa 02 
1,1205 


LEN OI 


zeroize coef 


delay 


left justify coef and adj exp 


examine sign of ans 


complement ens upper and lower if Sign is neg 


return to M.P. 
BEGIN 


zeroize temp storage 


test arguments as for. FMU EXCEPT DIVISOR =O 
gives error halt 
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1550 
1551 
1552 
1595 
1554 
1555 
1556 
1557 


"1560 
1561 
1562 
1563 
1564 
1565 
1566 
1567 


t950 
1512 
1572 
1212 
1574 
1212 
1576 
ОА 


1600 
1601 
1602 
160% 
1604 
1605 
1606 
1607 


Bono 
1611 
1612 
1615 
161} 
1615 
1616 
37 


1620 
1621 
` 1622 
1623 
1624 
1625 
1626 
1627 


0701 
6501 
0300 


6036 
2401 
6507 
5001 


2402 
4002 
0501 
2232 
6206 
2200 
0005 
0701 


6501 
0500 


6110 
4001 
2200 


Don 


0701 


6501 
0500 
1007 
0514 
4210 
2002 
3404 
6210 


5403 
6142 
0500 
Tite 
0000 
0000 
0000 
2200 


0002 
0701 
6501 
2002 
3404 
6305 
4002 
5112 


nzo OL 
nop 


lcn 14 
ste 10 
ldd 02 
spa Ol 
PIT 1O 


aod 03 
nzf +2 
1сп ОО 
hit 72 


ldc OO 


0002 
500 01 


n20 0L 


144 02 
sba 04 
njf 05 
std 02 
aob 12 


set counter 
coef check for start 


coef Ol greater than 02 jump to 1617 


jump to 1653 


delay 


jump to 1632 


in 1615 
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1630 
1631 
1632 
1655 
1634 
1635 
1636 
1637 


1640 
1641 
1642 
1643 
1644 
1645 
1646 
1647 


1650 
1651 
1652 
1653 
1654 
1655 
1656 
1557 


1660 
1661 
1662 
1663 
1664 
1665 
1666 
1667 


1670 
1671 
1672 
1673 
1674 
1675 
1676 
1677 


1700 
1701 
1/02 
2705 


oLol 
6206 
2200 
0002 
ОТОТ 
6501 
0300 
NES 


6220 
4724 
2002 
6310 
4402 
6721 
2200 
0004 


0701 
6501 
6613 
4402 
3404 
4002 
0401 
6630 


2200 
2001 
3403 
5001 
6752 
2350 
5002 
2551 


6107 
2200 
0005 
0701 
6501 
0500 
6005 
2401 


001 
2402 
1002 
7007 


Tann oL 
pjf 06 


14с 00 


0002 
500 Ol 
1120 01 
пор, 
aob 25 


pjí 2 
srb 2% 
laa 02 
1. бе 10) 
514 02 
njb 21 
ldc OO 
000% 


500 01 
nzb Ol 
рјо 13 
sra 02 
sbà Oh 
std 02 
dono 
pjb 30 


1а? ОО 
таа ол: 
sbd 03 
гад О1 
9619 2 
40:50 
std 02 
Job 51 


121 01 
14с OO 
0005 

son O1 
nzb Ol 
nop 

т 9 
led Ol 


stû Ol 
lcd. 02 
std 02 
jpi OT 


delay 


increase counter by one 


coef điv compl; jump to 1660 
on 1615 


jump to 1653 


jump to 1624 


jump to 1637 


form exponent 


place coef of ans in 0002 · 


check for sign of ans 


compl for neg ans 


return jump to М.Р. 





APPENDIX 11 


STATE GENERATOR AND TEST 


1. Theory. 

Linear approximation between samples, as follows: 

q+ SS А 11-1 

Both 9B and 8A are the average of four samples taken as quickly as 
possible (141.3 u sec/semple). The time, Qt, between 0B and OA is 
840 cycle times, or at 6.4 second/cycle time, 5,376 u sec. 
2. Memory requirements. 

Lower memory: 01-04, 07, 10, 20-47, and 61-67. 

Load State Generator: 2000-2322. 

Load Output Section: 2323-2357. 
3, Subroutines required. 

Constant Time Floating Point Arithmetic Package. (APP. I) 

A quick "Divide-by-4" routine, included with this appendix. 
4. The programs for the state generator and the "Divide-by-4" 


subroutine follow. Comments are included with the programs, 


А=23 





2010 


2011 
2012 
2013 
2014 
2015 
2016 
2017 
2020 


evel 
сога 
202% 
202% 
2025 
2026 


2950 


1500 
1101] 
7600 
4262 
1500 
1401 
1600 
4257 


7500 
1401 
1600 
4254 
1200 
1401 
1600 
4251 


4001 
0101 
Ообо! 
4007 
7061 
2001 
4021 
2240 


001 
0101 
0604 
4007 
uoo 
2001 
5021 
2227 


4001 
0101 
0604 
4007 
1061 
2001 
5021 
2216 


exf 00 


1401 

ina OO 
stf 62 
exf OO 
1401 

ina 00 
Ste 57 


exf 00 
1401 

ina 00 
зин 5 
exf 00 
1501 

ina 00 
STE Sl 


stá 01 
pta 


SE OF 


221 


37 
std 21 
laf hO 


std Ol 
pta 


SAO 


jp ol 
1G O1 
rad 21 
Дат 27 


std O1 
pta 


sta 07 
ірі 61 
хаа От 
гаа 21 
пара те 


STATE GENERATOR 





2031 
2032 
2055 
2034 
2055, 
2036 
2057 
2010 


2011, 
2052 
2043 
2044 
2045 
2056 
2017 
2050 


2051 
2052 
20902 
2054 
2055 
2056 
2057 
2060 


2061 
2062 
2063 
2064 
2065 
2066 
2067 
2070 


2071 
2072 
2075 
2074 
2919 
2076 
2077 
2100 


2101 
2102 
2105 
2104 
2185 
2106 
2101 
2110 


5001 


"0101 


0604 
5007 
7061 
2001 
5021 
2600 


0464 
0601 
6501 
0500 
6205 


OA2 
OAS 


OAL 

1200 
1401 
7600 
4256 
1500 
1401 
7600 


4253 
1900 
1401 
1600 
4250 
1200 
1401 
7600 


L245 
001 
0101 
060% 
4007 
1061 
2001 
4022 


2231 
4001 
0101 
0604 
НООТ 
1061 
2001 
5022 


sta 01 
pta 

adn Ol 
STA OT 


‚4р1 61 


lad 01 
rad 21 
Let 90 


làn 6h 


adn 01. 


nzb Ol 


‚пор 
pjf 0). 


0000 
0000 
0000 


0000 
exf 00 
1401 


апа 00 


stf 56 
exf 00 
1401 

ina 00 


OS 
exf 00 
1401 
1па 00 
stf 50 
exf 00 
1401 

ina 00 


stf 45 


pta 

adn Ol 
sta 07 
Joi ol 
144 Ol 
std 22 


laf sh 
544 01 
pta 

адл Ol 
58:01:01 
ІШІ Ol 
0 7 
rad 2 


GA avg in 0021 





2111 
2112 
2113 
2114 
2115 
2116 
2117 
2120 


2121 
22 
25 
212% 
2125 
2126 
ele | 
2130 


2131 
2132 
2152 
2134 
2155 
2156 
2157 
2140 


2141 
21h 2 
2143 
2144 
2115 
2146 
2147 
2150 


2151 
2152 
2153 
2154 
2155 
2156 
21517 
2160 


2161 
2162 
2163 
2164 
2165 
2166 
2167 
2170 


x US 


2223 
4001 
0101 
0605 
4007 
1061 
2001 
5022 


2212 
5001 
0101 
0604 
4007 
7061 
2001 
5022 


6206 
6306 
ОВ1 
OB2 
ОВ» 


0300 
5121 


4024 
OTOT 


4007 
202% 
7066 
2001 
4030 


2902 
1051 
20l, 
4003 
2045 
hook 


0101 


0604 


1007 
7065 
2001 
1032 
2002 
4033 
7200 
1401 


141 23 
514 01 
pta 

adn Oh 
stá 07 
арі 61 
ad 01 
red 22 


ldf 12 
sta 01 
pta. 

adn باون‎ 
sta 07 
101251 
20104071 
rad 22 


pjf 06 
аера 015 
0000 , 
0000 
0000 
0000 
nop 
504 21 


sta 21 
pta 

adn 09 
514 01 
194 24 
ірі 66 
144 01 
std 50 


adn Oh 


std OT 
jpi 65 
lad Ol 
std 32 
144 02 
std 33 
exf 00 
1401 


ӨВ аур іп 0022 


ӨВ - СА 


ӨВ - 04 to floating point 


OB DOT formed б 





2171 
2372 
2173 
2174 
2175 
2176 
2177 
2200 


0+00 


2202 
2203 
2204 
2205 
2206 
2207 
2210 


Bor) 
2212 
821% 
2211 
2215 
2216 
Е 
2220 


222] 
2222 
2225 
2221 


2226 


Ecc 
2250 


“2251 
2232 
садр 
2234 
2255 
2236 
Bor 
22h0 


2241 
2242 
2243 
2244 
2215 
2246 
2247 
2250 


1600 
1256 
1200 
1401 
7600 
4253 
1200 
1401 


7600 
4250 
7500 
1401 
7600 
4245 
001 
0101 


0604 
4007 
7061 
2001 
4025 
2234 
4001 
0101 


0605 
1007 
7061 
2001 
5025 
2223 
4001 
0101 


0608 
007 
1061 


2025 
2212 
5001 
0101 


0601 
007 
7061 
2001 
202% 


6506 
D 


ina OO 
stf 56 


exf 00 


1401 
ina 00 


5117:527 
exf OO ` 


1401 


ina 00 
stf 50 
exf OO 
1401 

ina 00 
stf 45 
std 01 
pta 


adn Ooh 


sta 07 


jpi 61 


302251 
184 01 
red 23 
931.06 
njf 06 
0000 


ОС аур іп 0025 





2251 
2252 
2529 
2254 
222 
2256 
2221 
2260 


2261 
2262 
2265 
2264 
2265 
2266 
2267 
2270 


2271 
2272 
2275 
2274 
2219 
2276 
221 
2300 


2301 
2302 
2505 


оса 
ОС5 
och 
0300 
3h22 
4025 
0101 
0605 


4007 
2025 
1066 
2001 
4034 
2002 
1035 
2044 


4003 
2045 
4004 
0101 
060% 
4007 
1065 
2001 


4036 
2002 
1037 


0000 
0000 
0000 
пор: 
604 22 
std 29 
pta 
adn 09 


ста 07 
144 25 
ірі 66 
144. 01 
std 34 
lad 02 
std 35 
iaa hà 


514 03 
144 45 


std Ol 


pta 

adn Oh 
sta 07 
ا‎ 
194 01 


std 36 
144 02 
Sa f 


"ес - OB formed 


ӨС - OB to floating point 


ӨС DOT formed 


continue to OUTPUT ROUTINE or CONTROL PROGRAM. 





SUBROUTINE DIVIDE BY FOUR 


1704 2001 144 01 pick-up argunent 
1705 6212 pjf 12+. "sign check 

1706 0 150 6 

1707 110 18 5 divide by با‎ 

1710 0102 is 1 

ШИГ 1200 ipe eo г 
+02۳ 

[1 5 5200 یپ‎ ۵ extend negative sign bit 
1714 6000 6000 

Ши: LOOL SG O store ansırar 

0001 9 l2 jump bz positive veth 
00۱ 15 © 

ШӘ) 0110 15 5 divide by L 

1721 0102 1s i 

1722 1200 lp2 00 f mask 

iow litt) lilt 

ies (0017 std 01 store answer 

1725 0300 nop 

1726 0300 nop time delay 

1727 0300 nop ۱ 

ШО ТООТ gpi OT return Jup ٦٣ 


Jump into and return jump out of this subroutine in 
the same manner as for the floating point routines. 





255% 
2351, 
2259 
2356 
5991 


2360 
2361 
2362 
2363 
2364 
2365 
2366 
2367 


2370 
2371 
2372 
2212 
2374 
2519 
, 2376 
2517 


200 
2h01 
2402 
2h05 
alkol 
2405 
2406 
2407 


2056 
001 
2057 
4002 
2200 


2002 
005 
2200 
6221 
LOOL 
0101 
0604 
007 


7065 
0101 
0604 
4007 
1067 
0500 


7500 


2401 
7304 
0021 
6103 
6002 
0020 
7010 
nom 


OUTPUT ROUTINE FOR e DOT 


BA ° 0 


144 56 


eco Or 
Tadaa, 
sta ве 
1ас ОО 


2002 
std 03 
пас со 
6221 
std 
pta 
adn 
514 


oh 


oh 
07 


Jpi 
pta 
adn 
std 
ps 
nop 
stä 
exf 


2.01 
out Oh 
0021 
8 O) 
2و‎ 2 
0020 
ОТО 
150201 


3۷۳162-۰61 پ7۷‎ ٣٢ 


6 DOT to fixed point 


A-30 




















| | D] 





